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Abstract. We present and analyze a novel concept for blue-detuned evanescent field 
surface traps for cold neutral atoms based on two-mode interference in ultra-thin 
optical fibres. When two or more transverse modes with the same frequency co- 
propagate in the fibre, their different phase velocities cause a stationary interference 
pattern to establish. Intensity minima of the evanescent field at any distance from the 
fibre surface can be created and an array of optical microtraps can thus be obtained 
in the evanescent field. We discuss three possible combinations of the lowest order 
modes, yielding traps at one to two hundred nanometres from the fibre surface which, 
using a few ten milliwatts of trapping laser power, have a depth on the order of 
1 mK for caesium atoms and a trapping lifetime exceeding 100 seconds. The resulting 
trapping geometry is of particular interest because atoms in such microtrap arrays will 
be coupled to any additional field propagating in the fibre via the evanescent field, 
thereby realising ensembles of fibre-coupled atoms. 
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Introduction 

Recently, the production of ultra-thin optical fibres with diameters smaller than the 
wavelength of the guided light has become possible in a number of laboratories [U EJ El H] . 
Such fibres have attracted considerable interest in the field of quantum optics due to their 
high potential for efficiently coupling light and matter [5J [3j E]. The guided modes in 
such ultra-thin optical fibres exhibit a unique combination of strong transverse confine- 
ment and pronounced evanescent field [7j. Furthermore, the strong radial confinement 
is maintained over the full length of the fibre waist, exceeding the Rayleigh range of a 
comparably focused freely propagating Gaussian beam by several orders of magnitude. 
This has been used in a number of experiments to couple atoms and molecules to the 
fibre mode via the evanescent field, showing that tapered optical fibres (TOFs) are a 
powerful tool for their detection, investigation, and manipulation: Recently, evanescent 
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field spectroscopy of a very small number of cold caesium atoms around a 500-nm diam- 
eter TOF has been performed [6]. In a similar experiment the fluorescence of resonantly 
irradiated atoms around a 400-nm diameter TOF, coupled into the fibre mode, has 
been detected and spectrally analyzed [31 [9]. The absorbance of organic dye molecules, 
deposited on a subwavelength-diameter TOF, has also been spectroscopically charac- 
terized via the fibre transmission with unprecedented sensitivity [10]. Finally, it has 
also been proposed to trap atoms around ultra-thin fibres using the optical dipole force 
exerted by the evanescent field [HI [12]. In this case, the atoms can be coupled to and 
trapped near a dielectric nanostructure without the need of additional external light 
fields. 

Here, we present a novel type of blue-detuned evanescent field trap for cold neutral 
atoms based on two-mode interference in such ultra-thin optical fibres. We consider 
a field-fibre configuration where only the four lowest order modes propagate, the fun- 
damental mode, HEn, and the first three higher order modes, TEoi, TMoi and HE21. 
If the modes are coherently excited, they will yield a stationary interference pattern 
while co-propagating in the fibre due to their different phase velocities. Previously, a 
similar scheme has been proposed for trapping atoms in the evanescent field of a two- 
dimensional planar waveguide. In this case, however, the interference of at least four 
waveguide modes is required in order to achieve three dimensional confinement of the 
atoms [8]. 

We explore all possible pairs of modes that can be used to trap cold neutral atoms 
in the intensity minima formed at the positions of destructive interference in the evanes- 
cent field surrounding the fibre. These combinations are HEu+TEqi, HE11+HE21 and 
TE01+HE21I3 We assume a cylindrical ultra-thin silica fibre with 400 nm radius for all 
three trapping configurations. The wavelength and the total power of the guided light 
as well as the power distribution between the modes have been chosen to fulfill the fol- 
lowing criteria for each individual trapping configuration: A three-dimensional trapping 
potential for caesium atoms, a depth of the trap on the order of 1 mK, and a trapping 
lifetime exceeding 100 seconds for an atom with an initial kinetic energy corresponding 
to a temperature of 100 /iK. 

This paper is organised as follows: Section CD is devoted to the analysis of the modal 
dispersion in ultra-thin optical fibres and to the presentation of the electric field equa- 
tions in the evanescent field for the three considered modes. In Sect. [2] the trap arising 
from the interference between the HE n and the TE i mode is presented (HEn-|-TE 01 ). 
Sections [3] and H] then treat the HE n -|-HE2i and TE 01 +HE2i traps, respectively. 

I The TM01 mode cannot be used to create a trap because its polarisation at the position of destructive 
interference cannot be matched with any of the other modes due to its large component in the direction 
of propagation. 
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1. Electric field and mode propagation in ultrathin optical fibres 

1.1. Mode propagation 

We consider a step-index optical fibre consisting of a cylindrical bulk of dielectric ma- 
terial with radius a and refractive index n±, surrounded by a second dielectric medium 
with infinite radius and refractive index n 2 . For the guided light in such a fibre a dis- 
crete set of propagation modes exists whose axial propagation constant (3 is fixed by the 
boundary conditions [13]. The number of modes and their axial propagation constant 
is then determined by the radius of the fibre a, the refractive indices of the two media, 
n\ and n 2 , and the wavelength of the light A via the parameter V = {2Tta/\)^n\ — nf. 

Figure [T] shows the axial propagation constant j3 for the first seven modes in the 
fibre normalized to the wavenumber in free space ko as a function of the V parameter. 
Note that /3/k lies between n 2 and n\ which is a condition that must be fulfilled by 
any lossless mode [13]. In the following, (3 will be referred to as propagation constant 
or phase velocity. The dashed vertical line located at V = 3.11 corresponds to the 
three configurations considered in this paper: An ultra-thin optical fibre of pure silica 
(ni = 1,452) with a radius a = 400 nm, surrounded by vacuum (n 2 = 1), and three 
similar wavelengths A = 849.0 nm, 850.5 nm and 851.0 nm. In this case, only four 
modes are allowed to propagate, the fundamental mode HE n and the first three non- 
fundamental modes, i.e., TE i, TM 01 , and HE 2 i. At this value of V the phase velocities 
of all modes differ significantly. This difference will cause an interference pattern to 
establish along the fibre and, in addition, results in different radial decay lengths of the 
evanescent field outside the fibre for each mode. The modal dispersion can therefore be 
used to create a tailored evanescent field resulting from the interference of two or more 
co-propagating modes. 



1.2. The fundamental HEu mode with quasi-linear polarisation 

The E field equations of the fundamental HEn mode with quasi-linear polarisation out- 
side the fibre, i.e., for r > a are given by [7]: 

E x (r,(f),z,t) = A n ^ n ^, a \ [(l - s n )K (q n r) cos(^ ) + 
2q u K^qna) 

+ (1 + Sii)K 2 {qxxr) cos(20 - <p )] exp[i(ut - (3 n z)] (1) 

Ey(r,<j),z,t) = A n ^ n ^ ^" [(l - s n )K Q (q n r) sm(ip ) + 
2<?n Ai(gna) 

+ (1 + s u )K 2 (q u r) sin(20 - y> )] exp[i{ut - (3 n z)] (2) 

E z (r,(f),z,t) = iA u J ^\ ll<X \ Ki(qiir) cos(0 - <p ) exp[i(ut - finz)] (3) 
K 1 {q u a) 
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where, 



Sll 



+ 



J[(h n a) K[(q n a) 



(h n a) 2 (g n a) 2 J L/i 11 aJ 1 (/i 11 a) quaK^qua) 



(4) 



hi = V^i-A 2 ! (5) 
q u = \J f5\ x - k\n\ (6) 

In the equations above, J'{x) (K'(x)) designates dJ(x)/dx (dK(x)/dx), a denotes the 
radius of the fibre, j3u the propagation constant of the HEn mode, and the angle <p 
gives the polarisation direction of the transverse electric field E± = (E x ,E y ), with 
fo = leading to x-polarisation and <po = n/2 to ^-polarisation. An is a normalisation 
constant of the fields that links the total power to the maximal field amplitude [15J. 
The quantity qu in Eq. (jSJ) is particularly relevant since it fixes the scale of the decay 
length of the fields outside the fibre which can be defined as An = 1/qn- The HEn is a 
hybrid mode since it is neither TE (transversal electric) nor TM (transversal magnetic) 
because the axial field components E z and H z are not zero [13]. The designation of 
quasi-linear polarisation stems from the fact that E z has a tt/2 dephasing with respect 
to E±, which results in elliptical polarisation except where E z = 0. 

Figure [2] shows a vectorial plot of the electric field component transverse to the 
fibre axis E± = (E x ,E y ) at t = and z = 0, with ipo set to zero. Note that the 
equations of the electric field inside the fibre used in this figure are not explicitly given 
here but can be found in [7]. The surface of the fibre is indicated by a gray circle. The 
calculations have been performed for a wavelength of A = 850 nm. The decay length of 
the evanescent field for these parameters is An = 164 nm. 

1.3. The TEqi mode 

We now present the E field equations of the TE i mode for r > a [13]. As can be seen 
from Fig. [TJ the TE i has, like the TM i, a cutoff value of V = 2.405. This is the lowest 
cutoff value of any non-fundamental mode and thus sets the single mode condition of a 
step-index optical fibre. 

Et(r,<t>,z,t) = ^ 'M^\ Boi K 1 (q 01 r) exp[i(cot - (3 m z)] (7) 



E z (r, 0, z, t) = E r (r, 0, z,t) = 



where, 



'01 



kln\ - fa (8) 



?oi = y^oi " (9) 
floi denotes the propagation constant of the TEoi mode and I?oi is the normalisation 
constant of the field amplitude. The TE 01 is classified as a transversal electric mode 
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since it has a vanishing ^-component of the electric field. Furthermore, the only non- 
vanishing electric field component is E^. The TE i thus possesses only one linearly 
independent polarisation state. The corresponding orthogonal polarisation state is the 
fibre eigenmode TM 01 . These two modes split due to the distinct influence of the fibre- 
vacuum boundary on the different E polarisation directions. 

Figure [3] shows a vectorial plot of the E field at t — and z = 0. The equations of 
the TE i mode inside the fibre used for this figure can be found in several text books, see 
for example [T3l IT4"] . The calculations have been performed for the same parameters as 
in Fig. As shown in Fig. [31 E vanishes at r = 0, which, together with the azimuthal 
symmetry of E^, produces a toroidal shape of the field amplitude distribution. The 
decay length of the evanescent field for the given parameters is Aqi = 277 nm. 



l.Jf.. The HE21 mode with quasi-linear polarisation 

The field equations of the fundamental HE 2 i mode with quasi-linear polarisation for 
r > a are given by [1 3j : 



where, 



E x (r, 0, z, t) = -A 21 p^ J ^ 2ia) [(l - 2s 21 )K 1 (q 21 r) cos(0 + 20 o ) + 
2q 2l K 2 {q 2 ia) 

+ (1 + 2s 2 i)K 3 (q 21 r) cos(30 + 20 o )] exp[i(ut - f3 2l z)} (10) 



E y (r,(j),z,t) = A 



21 



■[(1 - 2s 21 )K l (q 21 r) sin(</» + 20 o ) - 



E z (r,(p,z,t) = -%A 



2q 21 K 2 (q 21 a) 
1 + 2s 2 i)K z (q 2 ir) sin(30 + 20 o )] exp[i(u;t - (5 2X z) 

J 2 (h 21 a) 



21 



^2(^210 



-K 2 (q 21 r) cos(2(0 + O )) exp[i(ut - f3 21 z)\ 



S21 



+ 



J 2 (h 21 a) 



+ 



K' 2 (q 21 a) 



-(h 21 a) 2 (q 2 ia) 2 \ \-h 2l aJ 2 {h 21 a) q2iaK 2 (q 2 ia 



h 2 i = \ K n \ ~ P21 



921 



4$: 



k 2 n 2 



(12) 

(13) 

(14) 
(15) 



(3 2 i denotes the propagation constant of the HE21 mode and A21 is the normalisation 
constant of the field amplitude. 0o determines the polarisation direction of E±, with 
(f) = and 0o = tt/4 leading to two orthogonal polarisation states of the transverse 
electric field. The HE21 is a hybrid mode with six non-vanishing components of the E 
and H fields. Like for the HE 11; the designation of quasi-linear polarisation stems from 
the fact that E z has a tt/2 dephasing with respect to E±. 
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Figure H] shows a vectorial plot of the electric field component transversal to the 
fibre axis E± = (E x , E y ) at t — and z = with <p set to zero. Again, the equa- 
tions of the HE21 mode inside the fibre used for this figure can be found in several text 
books [T3J H3J- The calculations have been performed for the same parameters as in 
Fig. [2] with a polarisation direction given by 0o = 0. The decay length of the evanescent 
field for these parameters is A 2 i = 420 nm. 

2. HEn+TEoi trap 

We now show that an evanescent field surface trap for cold atoms can be obtained from 
the interference between the fundamental HE n mode, which we assume to be quasi- 
linearly polarised, and the TE i mode. By choosing the appropriate power distribution 
between the modes, an array of local minima of the field intensity at any distance from 
the fibre surface can be created at the positions where the two fields optimally can- 
cel. For blue-detuned light with respect to the atomic transition frequency a dipole 
force proportional to the negative gradient of the field intensity is then exerted on the 
atoms p2] confining them in the intensity minima. 

As an example, we discuss the properties of the above trap for caesium atoms. The 
trap can be created with 50 mW of light at a wavelength of 850.5 nm and the same fibre 
parameters as in Sect. HJ This power is realistic for such a fibre in vacuum: We could 
experimentally show that appropriately produced fibres with an even smaller radius of 
250 nm carry more than 300 mW of power in such conditions without fusing. With 72% 
of the power propagating in the HEn mode and 28% in the TE01 mode, a trap for cold 
caesium atoms with a trapping minimum at 134 nm from the fibre surface is formed. 
The depth of the trap is 0.92 mK and the trapping lifetime resulting from heating due 
to spontaneous scattering of photons exceeds 100 seconds for caesium atoms with an 
initial kinetic energy corresponding to 100 /iK. 

Figure [5] shows a contour plot of the trapping potential including the van der Waals 
surface potential [16] in the plane z = 4.61 /im. For the calculations, we use the van der 
Waals potential of an infinite planar silica surface [12]. The fibre surface is indicated 
by a gray circle and the equipotential lines are labelled in mK. The trapping minimum 
is located at <p = tc/2, r = 534 nm and z = 4.61 jum. The trapping minimum lies on 
the |/-axis because here the polarisation of the two modes matches and the interference 
is maximally destructive. This polarisation matching between the two modes can be 
understood when comparing Figs. [21 and [3j Note that while destructive interference 
takes place at = tc/2, there is constructive interference at <fi = 3tt/2. When varying 
the ipo parameter in Eqs. ([I]), (T5]) and ([3]), the polarisation direction of the HEn mode 
can be turned and thereby the azimuthal position of the trap can be varied because the 
potential has a cos(0 — fo) dependence. Using a harmonic potential approximation, we 
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Figure 1: Normalised propagation constant 
P/ko versus V parameter for the first seven 
modes in the fibre. The dashed vertical line is 
located at V = 3.11 which corresponds to the 
three trapping configurations considered in this 
paper. 
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Figure 2: Field plot of the electric field 
component perpendicular to the fibre axis E± = 
(E x ,E y ) for the HEn mode at t = 0, z = and 
for ip = (see Eqs. © and P)). The fibre 
is indicated by the gray circle. The following 
parameters have been used: a = 400 nm, n\ = 
1.452, n 2 = 1, and A = 850 nm. 
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Figure 3: Field plot of the electric field E for 
the TE01 mode at t = and z = (see Eq. ([J])). 
The fibre is indicated by the gray circle. The 
fibre parameters are identical to Fig. [2j 
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Figure 4: Field plot of the electric field 
component perpendicular to the fibre axis E± = 
(E x ,E y ) for the HE21 mode at t = 0, z = 
and for 4> = (see Eqs. (JTOJ) and ([II])). The 
fibre is indicated by the gray circle. The fibre 
parameters are identical to Fig. [2 
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calculate the azimuthal oscillation frequency to be lo^jlv: rs 1.07 MHz. The extension 
of the trap volume in the azimuthal direction for caesium atoms with a kinetic energy 
corresponding to 100 /iK is 34 nm. 

Figure [6] shows the contour plot of the trapping potential in the plane x — 0. The 
fibre surface is indicated by two vertical gray lines. The interference between the modes 
creates an array of traps in the axial direction with a periodicity given by the beat length 
of the two co-propagating modes, z = 27r/(/3n — /?oi) = 4.61 jum. In addition, there is 
a second array of traps on the opposite side of the fibre with same periodicity which is 
shifted by Zq/2. The potential has a sin((/3n — Poi)z) dependence in the axial direction. 
The axial trapping frequency is calculated to be u) z /2ir « 528 kHz. The extension of 
the trap volume in this direction for caesium atoms with a kinetic energy corresponding 
to 100 fiK is 68 nm. 

Figure [7] shows the trapping potential along the y-axis. The fibre surface is in- 
dicated by a vertical gray line. The solid black line corresponds to the sum of the 
light-induced potential and the van der Waals potential when 72% of the power propa- 
gates in the HE n mode. The dashed and dotted lines correspond to the same potential 
assuming slightly different power distributions between the modes: We define the pa- 
rameter t such that Pn = rP and Poi = (1 — r)P, where P denotes the total power 
transmitted through the fibre, Pn the power propagating in the HE n mode, and P i 
the power propagating in the TEoi mode. We assume that r can be controlled with a 
precision of a = 0.05a/t (1 — r ), i.e., a = 0.025 for r = 0.5. For the case of r = 0.72 
the power distribution between the modes r would then be controlled within ±0.022. 
We consider this value to be a conservative assumption for the precision of the power 
distribution between the two modes. For the case of Pn = (r + cr)P (dotted line) the 
trap is 27% shallower compared to the trap for Pn = tqP (solid line), whereas for the 
case of Pn = (t — a)P (dashed line) the trap is 30% deeper. While the trap depth 
increases when decreasing r, the trapping minimum is also shifted towards the fibre. 
When further decreasing r the depth of the trap thus drastically reduces because the 
van der Waals potential becomes larger than the light-induced potential. Furthermore, 
the potential barrier in the direction towards the fibre becomes narrower which would 
eventually lead to tunnelling of the atoms. The parameters presented here have been 
chosen in such a way that even with realistic experimental uncertainties the trap remains 
sufficiently deep and the tunnelling is negligible compared to the trapping lifetime. Note 
that the total potential is negative at its minimum due to the influence of the van der 
Waals potential. Since the ^-component of the electric field in the HEn mode vanishes 
at 4> — 7r/2, the polarisation in the two modes perfectly matches at the intensity mini- 
mum and the van der Waals potential at this position is the only influence on the atoms. 
We calculate the radial trapping frequency to be u r /2ir « 770 kHz and the extension 
of the trapping volume in the radial direction for caesium atoms with a kinetic energy 
corresponding to 100 fiK is 47 nm. 
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The calculations of the lifetime have been performed assuming caesium atoms with 
an initial kinetic energy equivalent to 100 fiK trapped in a three dimensional classical 
harmonic potential with oscillation amplitudes corresponding to the above given exten- 
sions of the trapping volume. Note that the trap is not perfectly symmetric along the 
y-axis (see Fig. [7j), we account for this fact by biasing the oscillation amplitude in this 
direction. From that, the mean squared field amplitude at the position of the atom 
has been calculated by integrating over all possible classical oscillation modes. Using 
this method, we find a scattering rate of 39 photons/second and a trapping lifetime of 
108 seconds. 



3. HEn+HE 2 i trap 

We now consider the trap arising from the interference between the HE n and the HE 2 i 
mode. It is created using 25 mW of light at a wavelength of 849.0 nm and the same fibre 
parameters as in Sect. [H The polarisation orientation of the modes has been chosen 
such that the trap forms at = 0. This corresponds to y?o = <Po — in Eqs. (CD) to (J3]) 
and (fTOl) to (fT2j) . respectively. With 84% of the power propagating in the HEn mode, 
i.e., r = 0.84 (see Sect. [2]), a trap at 152 nm from the fibre surface is formed. The depth 
of the trap is 1.2 niK and the trapping lifetime resulting from spontaneous scattering 
of photons exceeds 100 seconds for caesium atoms with an initial kinetic energy corre- 
sponding to 100 fiK. 

Figure M shows a contour plot of the trapping potential in the plane z = 3.45 /xm. 
The two dashed lines with their origin at the center of the trap indicate the two direc- 
tions with minimal potential barrier which, by consequence, determine the depth of the 
trap (see Fig. [TTj) . The trapping minimum is at = 0, r = 552 nm and z = 3.45 /im. It 
lies on the x-axis because here the polarisation of the two modes matches and the in- 
terference is maximally destructive. This polarisation matching between the two modes 
becomes apparent when comparing Figs. [2] and HI However, unlike the HEn+TE i trap 
considered in Sect. [21 the polarisation matching between the two modes is not perfect. 
This is due to the fact that the ratio E Z /\E±\ at the trapping minimum is different for the 
two modes and, therefore, the electric fields never cancel completely. Indeed, this stems 
from the orientation of E± at the position of the trap: When the transverse electric field 
is perpendicular to the fibre surface, a non-vanishing ^-component of the electric field 
arises [14]. This polarisation configuration results in a more intense evanescent field 
allowing the creation of a trap comparable to the one presented in Sect. [2] with only 
50% of the power. As a drawback, the intensity at the trapping minimum is not zero. 
When varying the parameters y?o and 0o in Eqs. ([TJ to (j3J) and (flDl to (1121) . respectively, 
the polarisation direction of the two modes can be rotated and thereby the azimuthal 
position of the trap can be varied. We calculate the azimuthal oscillation frequency to 
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Figure 5: Contour plot of the HEn+TE i trap in 
the plane z = 4.61 urn for the following parameters: 
P = 50 mW, r = 0.72, A = 850.5 nm, a = 400 nm, 
n\ = 1.452, and n% = 1. The fibre surface is indicated 
by the gray circle and the equipotential lines are 
labelled in mK. 
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be uJ<f,/27c pa 330 kHz. The extension of the trapping volume in the azimuthal direction 
for caesium atoms with a kinetic energy corresponding to 100 /iK is 104 nm. 

Figure [9] shows the contour plot of the trapping potential in the plane y = 0. Like 
in the HEn+TEoi trap, the interference between the modes creates an axial array of 
traps with a periodicity given by the beat length of the two co-propagating modes, 
Zq = 2tt/(Pu — P21) — 3.45 /an. Again, there is a second array of traps at the oppo- 
site side of the fibre with same periodicity and shifted by zq/2. The potential has a 
sin((/?n — @2i)z) dependency in the axial direction plus the offset due to the unbalanced 
^-components of the electric fields of the two modes. We calculate the axial trapping 
frequency to be w z /2n pa 610 kHz. The extension of the trapping volume in the axial 
direction for caesium atoms with a kinetic energy corresponding to 100 /iK is 58 nm. 

Figure [10] shows the trapping potential versus the position along the x-axis. The 
solid black line shows the radial trap for Pn = tqP, the dashed line for Pn = (to — cr)P, 
and the dotted line for P n = (r + cr)P, with r = 0.84 and a = 0.018. P denotes the 
total power propagating through the fibre and Pn the power propagating in the HEn 
mode. Again, r is assumed to be controlled with a precision of a = 0.05a/t (1 — r ). 
The light-induced potential does not vanish at the minimum due to the mismatching 
in the polarisation between the two modes. This leads to a higher scattering rate of 
57 photons/second compared to the trap presented in Sect. [2j The radial trapping 
frequency is u r /2iT pa 970 kHz. The extension of the trapping volume in the radial 
direction for caesium atoms with a kinetic energy corresponding to 100 /iK is 37 nm. 
Note that the depth of the potential shown in Fig. [10] does not correspond to the depth 
of the trap because, as mentioned above, the direction with minimal potential barrier 
for the atoms is not radial. Figure [11] therefore shows the trapping potential against 
the position along the direction with minimal potential barrier. The solid, dashed and 
dotted lines have been calculated for the same values of r as in Fig. [TO] We define the 
direction with minimal potential barrier I as the straight line that connects the potential 
minimum in the trap with the lowest local potential maximum. Note that I depends 
on t and has, per definition, its origin at the trapping minimum. Hence, the three 
minima of the potential profiles shown in Fig. [IT] are located at I = 0. The trap depth 
is then found to be 1.2 mK. For the case of Pn = (r + cr)P (dotted line) the trap is 
33% shallower compared to the trap for Pn = tqP (solid line), whereas for the case of 
Pn = ( r o — cr )P (dashed line) the trap is 17% deeper. Finally, we calculate a trapping 
lifetime of 106 seconds for caesium atoms with an initial kinetic energy corresponding 
to 100 /iK. Again, the tunnelling through the potential barrier in the radial direction 
towards the fibre (see Fig. [TOl is negligible compared to the lifetime of the atoms in the 
trap. 
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Figure 8: Contour plot of the HEn+HE 2 i trap in 
the plane z = 3.45 //m for the following parameters: 
P = 25 mW, r = 0.84, A = 849.0 nm, a = 400 nm, 
ni = 1.452, and ri2 = 1. The fibre surface is indicated 
by the gray circle and the equipotential lines are 
labelled in mK. 
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Figure 9: Contour plot of the HEii+HE 2 i trap in the 
plane y = for the same parameters as in Fig. [HJ The 
fibre surface is indicated by the two vertical gray lines 
and the equipotential lines are labelled in mK. 
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Figure 10: Plot of the trapping potential versus the 
position along the x axis for Pn = tqP (solid line), 
P n = (r + a)P (dotted line) and Pn = (r - cx)P 
(dashed line). The parameters are the same as in 
Fig. [8l The fibre surface is indicated by the vertical 
gray line. 



Figure 11: Plot of the trapping potential versus 
the position along the direction of minimal potential 
barrier l(r) for Pn = ToP (solid line), Pn = (to+<t)P 
(dotted line) and Pn = (tq — er)P (dashed line). The 
parameters are the same as in Fig. [HJ 
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4. HE 21 +TE i trap 

Finally, we consider the trap arising from the interference between the TEoi and the 
HE 2 i mode. It can be created using 30 mW of light at a wavelength of 851.0 nm and 
the same fibre parameters as in the above sections. The polarisation orientation of 
the modes has been chosen such that the trap forms at = 3tt/A and at <fi = — 7r/4. 
Note that this trapping configuration has two trapping minima in the same z-plane, 
whereas in the traps discussed before there is only one trapping minimum per z-plane. 
The polarisation orientation corresponds to 4> = in Eqs. ffTUj) . ffTTl) and f|T2|) for the 
HE 2 i mode. With 68% of the power propagating in the TE i mode, i.e., r = 0.68 (see 
Sect. |2]), a trap for cold caesium atoms with its trapping minimum at 184 nm from the 
fibre surface is formed. The depth of the trap is 1.4 mK and, like in the above cases, the 
trapping lifetime resulting from spontaneous scattering of photons exceeds 100 seconds 
for caesium atoms with an initial kinetic energy corresponding to 100 /iK. 

Figure [12] shows a contour plot of the trap in the plane z = 13.67 /im. Here, 
the trapping minima are shown to be at <fi = 37r/4, r = 584 nm, z = 13.67 /im and 
at <p = — 7r/4, r = 584 nm, z = 13.67 /im because the polarisation in the two modes 
matches at these positions. This is apparent when comparing Figs. [3] and HI Like for the 
HEn+TEoi case, the polarisation matching between the two modes is perfect because 
the ^-component of the electric field in the HE 2 i mode vanishes at the position of the 
trap. We calculate the azimuthal oscillation frequency to be lo^jl-n « 2.60 MHz. The 
extension of the trapping volume in the azimuthal direction for caesium atoms with 
a kinetic energy corresponding to 100 /iK is 14 nm. This strong confinement in the 
azimuthal direction stems from the behaviour of the polarisation of the electric field in 
the two the modes at the position of the trap. When increasing 0, the polarisation of 
the HE 2 i mode rotates clockwise, whereas the polarisation of the TE i mode rotates 
anticlockwise. This produces a fast polarisation mismatching between the two fields 
when displacing the position along the azimuthal direction and thereby a steep increase 
of the potential. 

Figure [TBI shows the contour plot of the trap in the zd— plane, where d — (y—x)/ y/2. 
The interference between the modes creates four axial arrays of traps with a periodicity 
of Zq = 27r/(/?oi — P21) = 13.67 /im. The two trapping minima shown in Fig. [T2l show 
the azimuthal positions of one pair of arrays. The second pair is shifted with respect to 
the first one by (f> = 7r/2 and z = zq/2. We calculate the axial trapping frequency to be 
uj z /2tc « 204 kHz. The extension of the trapping volume in the axial direction for cae- 
sium atoms with a kinetic energy corresponding to 100 /iK is 174 nm. This elongation 
of the trap compared to the traps presented in Sects. [2] and [3] stems from the large beat 
length between the TE i and the HE 21 mode. 

Figure [H] shows the radial trapping potential in the above defined zd— plane. The 
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solid black line shows the radial trap for P 01 = r P, the dashed line for P i = (r — cr)P, 
and the dotted line for P i = (r + cr)P, with r = 0.68 and o = 0.023. Again, r 
is assumed to be controlled with a precision of a = O.OSv/roCL-ro). For the case 
of Pn = (to + cr)P the trap is 25% shallower compared to the trap for Pn = r P, 
whereas for the case of Pn = (r — cr)P the trap is 36% deeper. Despite the vanishing 
light-induced potential at the trapping minimum, the total potential does not become 
significantly negative because the influence of the van der Waals potential at this dis- 
tance from the fibre surface is negligible. We calculate the radial trapping frequency to 
be u r /2n rs 770 kHz. The extension of the trapping volume in the radial direction for 
caesium atoms with a kinetic energy corresponding to 100 fiK is 47 nm. Since the beat 
length between the TE i and the HE 21 mode is large compared to the beat length in the 
other two traps, one would expect the radial size of the trap to be large as well. However, 
the difference in the decay lengths A21 — A01 is not the only factor that influences the 
radial profile of the trap. It is also determined by the exact functional dependence of 
the evanescent field for the different modes which results in a similar radial confinement 
compared to the HEn+TE i and HE n +HE 2 i configurations. Finally, we calculate the 
scattering rate and the trapping lifetime for caesium atoms with an initial kinetic energy 
corresponding to 100 fiK to be 62 photons/second and 114 seconds, respectively. 

5. Conclusions 

We presented three blue-detuned surface traps for cold atoms based on two-mode 
interference between the lowest order modes in the evanescent field around a 400-nm 
radius optical fibre. The trapping potential confines the atoms in all three dimensions: 
Radially, thanks to the tailored evanescent field, axially, thanks to the difference in the 
phase velocity, and azimuthal, because the different modes have different polarisation 
distributions around the fibre axis. The three traps have a depth of the order of 1 mK 
and a trapping lifetime that exceeds 100 seconds for caesium atoms with an initial kinetic 
energy equivalent to 100 fiK. Moreover, we have shown that the three traps are robust 
against experimental uncertainties in the power distribution between the modes. Such 
an array of optical microtraps in the evanescent field surrounding an optical fibre in 
combination with the highly efficient coupling of the atoms to the fibre modes provides 
a very promising framework, e.g., for experiments of storing and retrieving light with 
atomic ensembles. Finally, the selective excitation of the fibre modes is experimentally 
conceivable: One of the most promising methods is the generation of Gauss-Laguerre 
modes in free space that match the modes in the fibre jTT] . The mapping between 
the Gauss-Laguerre modes in free space and the modes in an ultra-thin optical fibre is 
therefore currently under investigation in our group. 
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Figure 12: Contour plot of the TE i+HE 2 i trap in 
the plane z = 13.67 um for the following parameters: 
P = 30 mW, r = 0.68, A = 851.0 nm, a = 400 nm, 
n\ = 1.452, and n<i = 1. The fibre surface is indicated 
by the gray circle and the equipotential lines are 
labelled in mK. 
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Figure 13: Contour plot of the TE i+HE 2 i trap 
in the zd— plane, where d = (y — i)/v2 for the 
same parameters as in Fig. [12j The fibre surface 
is indicated by the two vertical gray lines and the 
equipotential lines are labelled in mK. 
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Figure 14: Plot of the trapping potential versus the position along the d = (y—x) / y/2 
axis for Pn = tqP (solid line), Pn = (to + cx)P (dotted line) and Pu = (tq — c)P 
(dashed line). The parameters are the same as in Fig. 1121 The fibre surface is indicated 
by the vertical gray line. 
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